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We propose a theory of longitudinal resistivity in the normal phase of quasi-one-dimensional 
organic superconductors near the quantum critical point where antiferromagnetism borders with 
superconductivity under pressure. The linearized semi-classical Boltzmann equation is solved nu¬ 
merically, fed in by the half-filling electronic umklapp scattering vertex as derived from one-loop 
renormalization group calculations for the quasi-one-dimensional electron gas model. The momen¬ 
tum and temperature dependence of umklapp scattering has an important impact on the behaviour 
of longitudinal resistivity in the the normal phase. Resistivity is found to be linear in tempera¬ 
ture around the quantum critical point at which spin-density-wave order joins superconductivity 
along the antinesting axis, to gradually evolve towards the Fermi liquid behaviour in the limit of 
weak superconductivity. A comparison is made between theory and experiments performed on the 
(TMTSF) 2 PF 6 member of the Bechgaard salt series under pressure. 
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I. INTRODUCTION 

The quasi-one-dimensional (quasi-ID) organic conduc¬ 
tors (TMTSF) 2 X, called the Bechgaard salts, are known 
to be among the first examples of correlated electron sys¬ 
tems in which superconductivity (SC) borders with an- 
tiferromagnetism or spin-density-wave order (SDW) in 
their phase diagranP®. Both orderings are brought in 
close proximity by tuning pressure which can be achieved 
either hydrostatically or by chemical means. This results 
in a characteristic SDW-SC sequence of instabilities that 
shares many common traits with other unconventional 
superconductors including members of cuprat^^, heavy- 
fermiorP^, and pnictid^^ series of compounds. This 
sequence is typically preceded by a metallic phase with 
unusual properties. For the Bechgaard salts, an anoma¬ 
lous metallic phase occurs nearby the critical pressure Pc, 
as a quantum critical point (QCP) connecting the order¬ 
ing temperature Tc for SC at its optimal value with the 
vanishing Tsdw for The temperature depen¬ 

dence of longitudinal resistivity near this point is found 
to depart from the behaviour for a Fermi liquid. A 
detailed analysis has revealed that resistivity is instead 
linear in temperature close to Pc, as commonly found 
near a QCpI^. It smoothly evolves toward the Fermi liq¬ 
uid prediction as pressure is tuned away from Pc, in ap¬ 
parent correspon dence with the gradual suppression of 
Tc under pressurd- ^^ l ^^ l 

Generally linked to a QCP with SDW is the presence 
of spin fluctuations in the metallic state. These fluctu¬ 
ations may give rise to an important source of electron- 
electron umklapp scattering. As a mechanism of mo¬ 
mentum dissipation, umklapp enters as a key det ermi- 
nant in the temperature dependence of resistivit> ^^^ * ^^ [ 
In systems like the Bechgaard salts, this scattering pro¬ 
cess is primarily the consequence of a weak dimerization 
of the organic stacks which imparts a half-filled char¬ 


acter to the electron band and allows the longitudinal 
transfer of charge carriers across the Fermi surfac^^. 
The existence of SDW fluctuations in the metallic phase 
of these materials has been amply borne out by NMR 
measurement^I^®^. The temperature dependence of the 
nuclear spin relaxation rate shows a pronounced 

enhancement of the Curie-Weiss form, at variance with 
the linear—T Korringa law expected for a Fermi liquid. 
The amplitude of enhancement is in precise correspon¬ 
dence with the anomalous resistivity and the size of Tc 
over the broad range of pressure where superconductivity 
is presenfP^^, showing an intimate connection between 
SDW and Cooper pairing. 

From a theoretical point of view, one-loop renor¬ 
malization group (RG) studies on the quasi-ID elec¬ 
tron gas model with umklapp scattering have previously 
shown how the SDW-SC sequence of instabilities can 
be reproduced by varying the transverse next-to-nearest- 
neighbour hopping term , which acts as an antin esting 
term that simulates the role of pressure in the modeP^^. 
Along the antinesting axis, a critical value can be de¬ 
fined as the analogue of the actual QCP at Pc- This is 
where Tsdw is suppressed and a d-wave SC (SCd) insta¬ 
bility appears at a maximum Tc. The SDW-SCd pattern 
is also found to be accompanied by SDW fluctuations 
over a wide temperature domain of the non ordered state. 
Above for instance, the calculated SDW susceptibility 
in the normal phase fits a Curie-Weiss form in agreement 
with the one extracted from NMR experiments. This en¬ 
hancement of SDW correlations is mainly ascribed to the 
positive respo nse of umklapp scattering to the growth of 
d-wave pairin^^SEH^ 

From the one-loop electron-electron vertex functions 
obtained for the electron gas, the quasi-particle sc atter- 
ing rate has been computed from the self-energ}^^^^. 
Near the rate develops a linear T dependence above 
Tc, to which a Fermi liquid T^ component subsequently 
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adds and grows as t'j_ distances from the QCP and Tc de¬ 
creases. This qualitatively agrees with the polynomial T- 
analysis of resistivity data in systems like (TMTSF) 2 PF 6 
and (TMTSF) 2 C 104 over the w hole pressure domain 
where superconductivity is founcPSEH, An alternative 
self-energy approach to resistivity data has been pro¬ 
posed from a hot-cold spots picture of the SDW nesting 
along quasi-ID Fermi surface in the conventional frame¬ 
work of the QCP where superconductivity is neglectecP^. 

In this work we intend to push the calculation of longi¬ 
tudinal resistivity a step beyond the self-energy approach 
by solving the semi-classical Boltzmann equation. This 
is achieved numerically with the aid of the RG method 
which is employed to compute the momentum and tem¬ 
perature dependent umklapp scattering vertex function 
entering the collision term of the Boltzmann equation. 
Our procedure is similar in the outline with the one re¬ 
cently proposed by Buhmann et al^ in the study of 
transport properties of the hole-doped cuprate supercon¬ 
ductors. 

The results put forward below show that resistivity 
develops a metallic linear-T dependence as antinesting 
approaches its critical value where SDW and SCd or¬ 
ders meet. This quantum critical behaviour in resistivity 
ensues from the anisotropic growth of half-filling umk¬ 
lapp scattering with lowering temperature, a consequence 
of reinforcement of commensurate SDW correlations by 
SCd pairing in the metallic phase, which persists down 
to Tc. Above the reinforcement undergoes a grad¬ 
ual decline that defines an extended region of quantum 
criticality where the temperature dependence of resistiv¬ 
ity can fit a power law dependence p(T) ^ with an 
exponent a < 2. The exponent approaches the Fermi 
liquid limit a ^ 2 at low enough temperature and for 
sufficiently large namely where Tc becomes small. 
Anisotropy developed by umklapp scattering is found to 
have a sizeable impact on the momentum dependence of 
the scattering rate extracted from resistivity. 

The linearized Boltzmann transport theory is intro¬ 
duced in Sec. II and Appendix A. In Sec. Ill and Ap¬ 
pendix B, we present the RG approach of the quasi-ID 
electron gas, from which the results for the phase dia¬ 
gram along with the temperature and momentum depen¬ 
dence of umklapp scattering vertex entering the Boltz¬ 
mann equation are given. In Sec. IV, numerical results 
for resistivity and their analysis in terms of renormalized 
umklapp vertex function are detailed. A comparison with 
existing experiments performed on the (TMTSF) 2 PF 6 
member of the Bechgaard salts is presented. We sum¬ 
marize and conclude this work in Sec. V. 


II. BOLTZMANN EQUATION 

The semicalssical Boltzmann equation describes the 
variation of the quasi-particle Fermi distribution function 
/ due to collisions and external forces. In the presence 
of a spatially uniform and static electric field S coming 
from a one-body electrostatic potential, the Boltzmann 


equation reduces to the expression 


dfjk) 

dt 


— e£ ■ 'S/hkf — 


'dm- 

dt 


(1) 


where e is the electric charge. The collision integral for 
an array of Np chains of length L takes the form 


' dfjky 

dt 


{LNp)-^ Y, w{k,k2\km) 

k2,ks,k4 


x{/(fc)/(fe2)[l-/(fe3)][l-/(fe4)] 
-[l-/(fc)][l-/(fe2)]/(fc3)/(fe4)}, (2) 


for scattering in and out processes. The electron-electron 
(longitudinal umklapp) contribution to scattering rate w 
is obtained from the Fermi golden rule, 

1 2 
w{k,k2;ks,k4) = -|(/c,/C2|^3|^3, ^ 4 ) - AJ2|^3|^4, A^s)! 

X ^k-\-k2,k3-\-k4-\-G^^{£k + ^^2 “ ^^3 ~ ^k^) 


(3) 


where G = (Akp^O) is the longitudinal reciprocal lat¬ 
tice vector for half-filling umklapp scattering, and kp = 
7r/(2a) is the longitudinal (ID) Fermi wave vector for 
dimerized chains. 

To linearize the Boltzmann equation, we introduce the 
normalized deviation function in the Fermi distribu¬ 
tion functiorP^, 


f{k) = 


1 


(4) 


where /3 = 1/T {kp = 1). The electron energy spectrum 
is the one of the quasi-ID electron gas model 


6^ =hvp{pk - kp) -h e±{kj_). (5) 


The spectrum comprises a longitudinal part linearized 
around the right/left ID Fermi points pkp = :Tkp with 
vp as the longitudinal Fermi velocity. The transverse 
part is given by 


^±{k_\_) = —2t± cos k_\_d± — 2d_^cos2k±d±, (6) 


where t± and are the first and second nearest- 
neighbour interchain hopping terms and dp is the in¬ 
terchain distance. 

The first order expansion f{k) f^{k) + f^{k)[l — 
/^(/c)]0fc of the distribution function yields after substi¬ 
tution in 0’ 

C(pk = eps ■ Vk, (7) 

where £ ^ 0{(j)). This is the linearized Boltzmann equa¬ 
tion in which the collision operator C obeys the relation 

£4>k = ( 8 ) 

k' 

The matrix elements are of the form 
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+ 5k2,k' ~ ^k3,k' — 5ki,k') — 


{LNpf k,MM‘^ 

/0(fe2)[l-/°(fc3)][l-/°(fc4)], 


[l-/°(fe)] 


i=l 


which separates into a diagonal and three off- 

diagonal terms In the Appendix A, each term 

is evaluated explicitly in the framework of the quasi-ID 
electron gas model introduced in Sec. |III[ 

For an electric field S = Sa^ oriented along the chain 
direction, the corresponding electrical current density 
ja = allows US to extract the conductivity ( 7 ^, as 

the inverse of longitudinal resistivity To first order 
in 0 , the current density is given by 

k 

- - f{k)]4)k^. (10) 

k 


Since for all temperatures of interest, the product /^[l — 
/^] is strongly peaked at the Fermi level, the deviation 
(j)h = (t)k±_ can be evaluated at the Fermi wave vec¬ 
tor whose location on the p Fermi 

surface sheet is entirely parametrized by k± from the 
equation e^p = 0. By introducing the normalized devi- 

K jp 


ation (l)kj_ = 4>k I / {P e£ad±) and after an energy integra¬ 
tion following (A4), the 2 D conductivity or the inverse 


resistivity in the longitudinal direction becomes 

o2 




-1 


= -^{4>k^)FS, 


( 11 ) 


where ( 0 /c^)fs = ^P^k± deviation func¬ 

tion averaged over the Fermi surface. The function 0/.^ 
satisfies the equation 


iik'j^ 


k± ,kj 


0/c' — I 5 


( 12 ) 


whose explicit expression is given in (A6) and for 


which we have introduced the dimensionless operator 
£[ 2 ] _ above expression is fed by the mo¬ 

mentum and temperature dependence of the umklapp 
vertex function. This function is provided by the RG 
approach to the quasi-ID electron gas model that is in¬ 
troduced next. 


III. RENORMALIZED UMKLAPP 
SCATTERING FOR THE ELECTRON GAS 
MODEL 


A. The model 


The quasi-ID electron gas model in its standard form 
comprises, besides its one-electron energy spectrum 


three electron-electron coupling constants that will be 
defined here on the Fermi surface sheets These 

are the backward and the forward scattering amplitudes 

^^, 4 ) 5 5 ^^, 4 ) 

normal scattering processes between right and left mov¬ 
ing carriers, and P3(A;^ 1 , A?^ 2 ^ 3 ’^^ 4 ) umklapp 

scattering along the chains, whose momentum conserva¬ 
tion involves the reciprocal lattice vector G = 0) at 

half-filling. These coupling constants can be seen as phe¬ 
nomenological parameters of the model whose range at 
the bare level can be fixed from experimental data. We 
will follovPSHUl^ 

and use the typical values (normalized 
by h'KVp)^ gi — 0.32 and g 2 — 0.64 for the two normal 
processes, consistent with the observed enhancement of 
uniform magne tic su sceptibility!^ and the scale of Tsdw 
at low pressur^^^m^ As for the umklapp term, its bare 
value is primarily proportional to the dime rizati on gap 
Aj^of the organic stacks (^3 « which 

is small in the Bechgaard salts. Although rather weak at 
the bare level, umklapp processes have a strong influence 
on spin fluctuations at low energy and enter as a key de¬ 
terminant in the temperature dependence of longitudinal 
resistivity as we will see. The bare umklapp amplitude 
will thus be varied in the interval g^ — 0.02... 0.035, in 
accordance with the range of values expected for AjSl. 
For the model spectrum we shall use Ep = 3000K 
and t± = 200K for the Fermi energy and transverse hop¬ 
ping, as representa tive fi gures of the band structure for 
the Bechgaard salt d^^ * ^^ *. The transverse second nearest- 
neighbour hopping which introduces nesting alter¬ 
ations, will be varied continuously to simulate pressure 
effects. 

We follow the lines of previous works and apply the 
Kadanoff-Wilson RG approach to the quasi-I D electron 
gas model which is outlined in Appendix 
Each constant energy surface is first divided into 
Np{= 60) patches, each centered on a discrete value of k± 
parametrizing each Fermi sheet = {k^{k±)^k±). The 
successive integration of electronic degrees of freedom as 
a function of the energy distance from the Fermi sur¬ 
face comprises perturbative contributions to the scatter¬ 
ing amplitudes coming from closed loops, vertex correc¬ 
tions and ladder diagrams of the density-wave (Peierls) 
and electron-electron (Gooper) scattering channels. The 
resulting RG flow equations for the -dependent scat¬ 
tering amplitudes ^ 1 , 2 , 3 (^) 4 !, k ± 2 ] k±s^ k± 4 ) on the Fermi 
surface are reproduced in ( |B 1 | ) of Appendix B. Their k± 
dependence takes into account the transverse momen¬ 
tum variables in both Peierls and Gooper channels, along 
with the longitudinal momentum transfer in each chan- 
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FIG. 1. (Color online) Calculated phase diagram of the 
quasi-one-dimensional electron gas model showing the SDW- 
SCd sequence of instabilities as a function of antinesting 
and initial gs. 


nel which is adjusted so that the momentum of ingoing 
and outgoing particles fall on the warped Fermi surface. 

These flow equations are co mpl eted by those of the sus¬ 
ceptibilities Xh{Qo) given in (B4). These serve to signal 
an instability against either /i =5DW or /i =SCd types of 
ordering at the corresponding wave vectors = ( 2 /ci?, tt) 
and Qq = 0 . 
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FIG. 2. Renormalized umklapp scattering amplitude pro¬ 
jected in the {k±i,k± 2 ) plane in the normal phase at dif¬ 
ferent temperatures for the initial gs = 0.03 {t'* = 38K). (a): 
4 = 33K (< 4); (b) : 4 = 44K (> 4). 


B. The phase diagram 

The integration of the flow equations for the scatter¬ 
ing amplitudes ( |B 1 [ ) and susceptibilities ( B4|B5 ) leads to 
their corresponding renormalized values at temperature 
T. For the above set of model parameters, the singulari¬ 
ties in XsDw(2/ci7,7r) and Xscd(O) at temperatures Tsdw 
and Tscd enable to follow the SDW-SCd sequence of in¬ 
stabilities as a function of the antinesting parameter 4 
and initial g^. The resulting phase diagrams are shown 
in Fig. The variation of Tsdw shows a characteristic 
monotonic decrease with up to the approach of the 
critical value 4 where Tsdw undergoes a rapid decrease, 
followed by the emergence of a SCd instability. The lat¬ 
ter reaches its maximum Tc at , followed by its steady 
decrease for larger The impact of increasing the ini¬ 
tial gs on this sequence is to magnify antiferromagnetism 
and Tsdw: which ultimately translates into an increase 
of the critical and an enhancement of d-wave Cooper 
pairing, as shown by the upward shift of the Tc line in 

Fig- [3 


C. Umklapp scattering 

In the SDW sector of the phase diagram, the umk¬ 
lapp scattering amplitude gs develops like xsdw a singu¬ 
lar growth at the approach of Tsdw- For and 

different temperatures. Fig. [^a shows the contour plot 


of g 3 {k±i^—k±i^k± 2 ) projected in the {k±i^k± 2 ) plane 
at zero transverse momentum for the pairs of ingoing 
and outgoing particles. As shown in the Figure [^a, 
the maximum of scattering intensity is mainly concen¬ 
trated along the lines k ±2 = k±i Ttt, which is congruent 
with the transverse component of the SDW correlation 
wave vector qq = ( 2 /ci 7 , 7 r). Along these lines, the sin¬ 
gular growth is the strongest around the spots ( 0 , Ttt) 
and (±7r,0), whose components differ from the points 
k± = ±7r/4,±37r/4 on the Fermi surface where nesting 
is predicted to be optimal for the spectrum (|5|. The 
fact that the warmest regions of scattering on the Fermi 
surface differ from the nesting prediction reflects the in¬ 
fluence of Cooper scattering channel in the anisotropic 
flow of ^3 in Eq. (Bl). 


If we now consider the low temperature variation of 
umklapp in the SCd sector of the phase diagram at 
we note according to Fig. that gs is no 
longer singular for all temperatures down to Tc. How¬ 
ever, as temperature is decreasing we still observe an 
enhancement of g^ along the SDW lines of scattering 
k ±2 = k±i ± TT. These are associated with peaks at 
k± = 0, ± 7 r on the Fermi surface whose amplitude in¬ 
creases with the bare gs^ as shown in Fig.[^ According to 
the combination ^sdw = 92 T ^3 of couplings entering in 
the SDW susceptibility [see Eq. (B 6 )], a positive increase 
in umklapp goes hand in hand with an increase of SDW 
spin correlations. This has been shown to yield a temper¬ 
ature dependence for Xsdw(Qo) ^ l/(T-l-0) of the Curie- 
Weiss form in the same temperature domain, with an en- 
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a) ^^3=0.02 , ti=35 K 



FIG. 3. (Color online) Renormalized umklapp scattering 
amplitude projected in the {k±i,k± 2 ) plane in the normal 
phase of the SCd sector of the phase diagram,as obtained at 
different temperatures and initial values of gs. 

ergy scale 0 for SDW correlations approa ching ze ro at 
and raising linearly with t'j_ from the Occur¬ 

ring despite poor nesting conditions, the enhancement is 
the consequence of a positive feedback of d-wave Cooper 
pairing on SDW correlations. As mentioned earlier, this 
is made possible by a -dependent coupling of gs to nor - 
mal scattering processes g 2 and gi. According to ( |B1[ ), 
the first Fourier component in cos(/cxi(i^) cos{k± 2 d±) of 
these couplings grows as temperature is lowered to finally 
become singular at the approach of Tc where Xscd(O) is 
divergent. 

Anisotropy of umklapp scattering in momentum space 
along with its temperature dependence enter as key fac¬ 
tors in the determination of resistivity at low tempera¬ 
ture. 


IV. NUMERICAL RESULTS FOR THE 
BOLTZMANN EQUATION 


We first consider in Fig. the general features of lon¬ 
gitudinal resistivity over the whole temperature interval. 
This includes the high - effectively ID - temperature do¬ 
main at T > t_L, its crossover to the 2D region at T < t± 
where the warping of the Fermi surface becomes quan¬ 
tum mechanically coherent, and finally the temperature 
domain < T < where antinesting effects come into 
play at the approach of instabilities. 

From the Figure we see that in the ID tempera¬ 
ture region of both SDW and SCd sectors of the phase 
diagram, the resistivity is not metallic but grows with de¬ 
creasing temperature. This is the result of relevant umk¬ 
lapp scattering whose flow in the ID domain is character¬ 
ized by essentially no k± dependence, that is, no trans¬ 
verse correlations, as illustrated in Fig. [^b (top panel). 
The growth completely overcomes the linear T-resistivity 
decrease obtained in the limit of a constant gs (continu¬ 
ous line of Fig.|^. This insulating behaviour is consistent 
with previous ID results obtained at half-fillin ^^^ l ^^ i 

When t'j_ < T < t_L, the system crossovers to its 2D 
regime with antinesting effects remaining weak. In this 
temperature domain, the relevance of g^ is still strong 
enough to outweigh the crossover toward the low tem¬ 
perature Fer mi liqu id In T-resistivity behaviour found 
at constant g^^^^ (continuous line of Fig. |^. Resistiv¬ 
ity remains insulating in this region. It is only when T is 
further lowered below that the anisotropic growth of 
^3 is significantly damped. At in the SCd sector 

where umklapp, though still increasing, is no longer sin¬ 
gular at finite T, resistivity goes through a maximum at 
T ^ and is followed at lower temperature by a metal¬ 
lic behaviour. Metallic behaviour is also possible above 
^SDW ill the SDW sector of the phase diagram provided 
is relatively close to the critical as portrayed in 
Fig. a As we will show next, however, metallicity is in 
general not equivalent to a Fermi liquid behaviour. De¬ 
viations from the latter are found and correlated to the 
distance from ; it is only at strong ^ that Fermi 
liquid features tend to be restored at sufficiently low tem¬ 
perature. 


B. Resistivity at quantum criticality and beyond 


We can now proceed to the numerical solution of the 
linearized semi-classical Boltzmann equation ([T^, using 
(A6) and the above RG results for the renormalized umk¬ 
lapp vertex part. The solutions for the normalized func¬ 
tions have been obtained for Np = 60 patches for 
the Fermi surface sheet k^{k±), which assures a conver- 
gen ce f or the solutions of the Boltzmann equation (12) 
or (A6) that is independent of the discretization of the 
transverse Brillouin zone. Following the tempera¬ 
ture dependence of longitudinal resistivity can then be 
obtained by tuning the values of t'j_ and initial g^. 


We now focus on the low temperature profile of pa at 
< T < on either side of the QCP Starting with 
the SDW side of the phase diagram, the Fig. [^displays 
the temperature dependence of resistivity in the normal 
phase of the SDW instability at different and bare gs 
values. When t'j_ is relatively far below resistivity 
remains insulating like with dpa/dT < 0 down to Tsdw- 
The growth of umklapp scattering and therefore of SDW 
fluctuations is sufficiently strong to maintain the source 
of inelastic scattering and the insulating response of the 
electron gas. In these conditions, the behaviour varies 
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FIG. 4. (Color online) Calculated longitudinal resistivity over 
the whole temperature interval for different 9s = 0.03 

{t'Z = 38K). The continuous line refers to the Fermi liquid 
limit at constant scale-independent gsi^ 0.03) yielding linear 
T-resistivity in the high temperature ID domain (T > t±) 
followed by a crossover toward a 2D In T regime at low 
temperature. 

weakly with the initial amplitude of 

The temperature profile of resistivity qualitatively 
change at the approach of At T ^ pa reaches a 
maximum followed by a crossover to a metallic behaviour 
which persists down to very close to Tsdw- The sin¬ 
gularity in umklapp, being still present despite sizeable 
nesting alterations, generates a resistivity upturn due to 
critical scattering, as meant in Figs and In the 
metallic regime, however, the resistivity deviates from 
the Fermi liquid prediction. A suitable way to describe 
these deviations is to express resistivity as a power law, 
pa{T) ^ with a local exponent a{T) defined by 

aiT) = ^\np,{T). (13) 

The values of the exponent a are inserted in the calcu¬ 
lated phase diagrams and shown in Fig. for different 
bare amplitudes of g^. Thus in approaching from be¬ 
low the exponent a{T) crossovers from negative (insulat¬ 
ing) to positive (metallic) values. Sufficiently close to 
Tsdw is small and a{T) reaches unity for T < lOK, which 
corresponds to the linear-T resistivity usually associated 
to the QCP. 

This linearity can be qualitatively understood if one 
considers that at low temperature resistivity is roughly 
speaking proportional to T^(^ 3 )fs- This would indicate 
that the average of the square of umklapp scattering over 
the Fermi surface, (^Dfs ^ 1/T, essentially grows as 
the inverse of temperature above the transition. Note 
that this temperature dependence is similar to the SDW 
susceptibi lity th at goes like xsdw ^l/TatT^|0|in 
this regioiP^^. 

If we now consider the SCd region at the 




FIG. 5. (Color in line) Calculated resistivity at low temper¬ 
ature in the SDW domain of the phase diagram for different 
tY and ^3- 


linearity continues to be seen over a finite region of the 
phase diagram. From Figures and it can be found 
over a large temperature interval above Tc near . How¬ 
ever, a shift in the value of ol{T) toward larger values 
takes place at low enough temperature especially if g^ is 
relatively small initially. As a function of temperature a 
fan shape like region then unfolds from the neighbour¬ 
hood of in which o:{T) 1. This region of linear 

resistivity in the phase diagram broadens as the bare 
^3 increases and spin fluctuations in the normal phase 
and then Tc are enhanced. Thus contrary to the situa¬ 
tion expected for a classical QCP tied to SDW ordering 
alone^, this region of linearity is not quickly followed as 
a function of t'j_ by a -resistivity. The crossover re¬ 
gion 1 < a{T) < 2 where deviations are present is rather 
wide and covers the whole region where superconductiv¬ 
ity is present, as portrayed in Figure for different gs. 
It is only in the low temperature domain close to Tc and 
for t'_^ well above that the boundary with the Fermi 
liquid limit spears and -resistivity becomes visible 
(see Figures]^ and [^, namely where (^ 3 )fs becomes es- 
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FIG. 6. (Color online) Calculated resistivity at low temperature in the SCd domain of the phase diagram for different and 
93- 


sentially temperature independent. Similar features were 
found in previous RG calculations of the electron-electron 
scatterin g r ate that included both normal and umklapp 
processed. These were attributed to the presence of the 
superconducting instability that still influence the growth 
of the coupling constants; it gives rise to a positive feed¬ 
back of SCd pairing on umklapp scattering which can 
account for this extended quantum critical behaviour. 


C. Anisotropy of the scattering rate 

It is instructive to extract from the above results for re¬ 
sistivity some information about the momentum resolved 
- /c^-dependent - scattering rate. According to ( pT] ), the 
conductivity <7^ = (JkA_i is a sum of contributions com¬ 
ing from each patch. One can associate to each contribu¬ 
tion a patch resistivity as an element of 

an array of parallel resistors proportional to a scattering 
rate at kFp{k^i) on the Fermi surface. 


The Figure]^ shows the variation of as a function 
of k^i. In the SCd regime, there are peaks at k± = 0 
and ±7r indicating that the scattering rate is larger at 
these points. These warmer regions of scattering are the 
same as those of umklapp scattering amplitude along the 
Fermi surface in Fig. and then differ from the positions 
k±_ = ±7r/4 and ±37r/4 attributed to the best nesting 
conditions for the spectrum. The resulting momentum 
profile of over the Fermi surface is also congruent 
with the one previously found by a direct RG calculation 
of the scattering rate from the one-particle self-energy at 
low temperature, which takes into account both normal 
and umklapp processe^^. As stressed in RefP these 
warmer regions of scattering turn out to coincide with 
the locations of the maximum for a d-wave gap along the 
Fermi surface. Since a single scattering channel (Peierls) 
calculation would invariably leads to the best nesting 
points as the warmest spots of scattering on the Fermi 
surface, it follows that one can ascribe the shift in the 
spots on the Fermi surface to the presence of the interfer- 
































FIG. 8. (Color online) Variation of resistivity components 
along the Fermi surface at different temperatures: SDW 
(top) and SCd (bottom) domains. 


FIG. 7. (Color online) Variation of the power law exponent 
Oi of resistivity in the phase diagram for gs = 0.026 (top) and 
gs = 0.033 (bottom). The arrows indicate the location of 


ing C ooper channel in the complete one-loop calculation 
(Bl). The interference between both scattering channels 
being present in the SDW sector of the phase diagram as 
well, the anisotropy of pkj^ is found to be similar above 
^SDW, as displayed in Fig. 


D. Comparison with experiments 

In the following we shall make some comparison of the 
above results with experiments. Starting with the SDW 
sector of the calculated phase diagram at low nesting de¬ 
viations, the results displayed in Figures and show 
the absence of metallic resistivity above Tsdw- This is 
at variance to what is actually observed in system like 
(TMTSF) 2 PF 6 series of the Bechgaard salts series be¬ 
low Pc where a nearly behaviour for resistivity is 
well known to be observed above the resistivity upturn 
at Tsdw - (see also Fig. 1^. The tem¬ 

perature profile obtained for resistivity well below is 
rather characteristic of the one seen in the cousin sulfur 


based compounds (TMTTF) 2 X (X=PF6, Br, ...) at low 
pressur^^H^ which are characterized by a ID Mott insulat¬ 
ing gap at high temperature due to umklapp scattering, 
followed at lower tempera ture by a transition to an anti¬ 
ferromagnetic ground stat cj^^FS JMJ^ Under moderate pres¬ 
sure, the normal phase of these systems becomes metallic 
with a power law behaviour down to Tsdw which would 
agree qualitatively with the calculated positive values for 
a found at low temperature not too far below in Fig- 
ure[^ However in contrast to the results shown in Fig.[^ 
no maximum in resistivity is found experimentally for 
longitudinal resistivity at a scale that would correspond 
to These discrepancies provide some indication 

that inelastic scattering at intermediate temperature and 
above would come from processes other than half-filling 
umklappP^®^. 

On the SC side, above the QCP at Pc, detailed resis¬ 
tivity measurements for (TMTSF) 2 PF 6 have been per¬ 
formed in the low temperature domai n up t o about 15 
times Tc {Tc of the order of 1 K at p^)| iQPd . The resis¬ 
tivity data are reproduced in Figure At the pressure 
of 11.8 kbar, close to the Pc — 9kbar, resistivity is found 
to be T—linear {a = 1) from up to lOTc down to Tc and 
even below when a small magnetic field is applied, as 
shown in the log-log plot of Figure [^b. As pressure is in¬ 
creased, the linearity decreases in strength and the resis- 
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tivity develops some upward curvature which translates 
into a gradual increase of the power law exponent a, but 
which still remains smaller than the Fermi liquid limit 
(a = 2. It is only at much larger pressure (P = 20.8 kbar 
in Fig.|^, when Tc becomes small, that the behaviour 
tends to be retrieved as shown in Fig. [^b. Deviations to 
the Fer mi liqu id behaviour were found to scale with the 
size of 

These features are consistent with the results of Fig.[^ 
where deviations from linearity develop with increasing 
the antinesting parameter so that the resistivity exponent 
reaches 2 when the superconducting phase is about to 
disappear. The calculated evolution of the exponent a 
near and above is in qualitative agreement with the 
results of Fig. 

The theory, however, fails to predict the restoration of 
an apparent Fermi liquid T^-behaviour in the interme¬ 
diate temperature region, namely above the domain of 
anomalous resistivity. As mentioned above, this Fermi 
liquid behaviour is also found in the SDW sector of the 
phase diagram above Tsdw ^ lOK (see Fig. |^. In the 
present approach, this would indicate that the system 
crossovers to a state where umklapp renormalization be¬ 
comes essentially absent, in spite of reinforced nesting 
properties resulting from either a temperature increase 
or a reduction of . These features are at variance with 
the results of the present model that invariably predicts 
the reverse. For systems like the Bechgaard salts with 
one type of carriers tied to a single component Fermi sur¬ 
face, these features may indicate that a ‘colder’ source of 
inelastic scatterinj^, not included in the present calcula¬ 
tions, gains in importance at higher temperature. It sur¬ 
passes around lOK the anomalous power-law behaviour 
resulting from the sole basis of electronic half-filled umk¬ 
lapp scattering. This is also inferred from the calculated 
values of 3D resistivity for the present model. The ac¬ 
tual resistivity is obtained from the product of the 2D 
part by the interchain lattice constant c, namely paC. 
For the Bechgaard salts, c = 13.52A1^ and the resistiv¬ 
ity can then be expressed in pCl • cm units, as shown in 
the right-hand side scale of Figures |4|[6l The comparison 
with the Fig.j^b for the inelastic part of resistivity indi¬ 
cates that the data reach one order of magnitude larger 
values at lOK, that is where the change of regime occurs 
(see Fig. [^. 


V. SUMMARY AND CONCLUDING REMARKS 

We have studied the temperature and pressure depen¬ 
dence of electrical resistivity of the Bechgaard salts quasi- 
one dimensional organic superconductors in the vicinity 
of their quantum critical point where superconductivity 
borders on spin-density-wave order under pressure. We 
have proceeded to the numerical solution of the semi- 
classical Boltzmann equation in which the collision umk¬ 
lapp electronic term has been computed by one-loop RG 




FIG. 9. (Color online) (a): Longitudinal resistivity of Ref.^ 
observed at different pressures or Tc for (TMTSF) 2 PF 6 ; (b): 
Inelastic part of the data (Apa = Pa — Pa) traced in a log- 
log plot showing the growth of the power law exponent a for 
resistivity with pressure. 

calculations in the framework of the quasi-ID electron 
gas model in the presence of tunable nesting alterations. 

The temperature dependence of longitudinal resistiv¬ 
ity pa has thus been obtained over the whole temperature 
domain from the ID high temperature region, which is 
invariably insulating, down to the different effective 2D 
regimes where either an insulating SDW or SCd instabil¬ 
ities are found depending on the strength of antinesting 
parameter Antinesting proves to be the characteris¬ 
tic scale of the model for weakening the growth of umk¬ 
lapp term under scaling, allowing the onset of anoma¬ 
lous metallic resistivity near the quantum critical point 
along the antinesting axis. Marked deviations from 
the Fermi liquid T^-behaviour have been found and de¬ 
scribed in terms of power law resistivity pa ^ in the 
chain direction. Linear resistivity is found within a char¬ 
acteristic fan shape region of the phase diagram near , 
followed by a gradual increase of the exponent a, reach¬ 
ing the Fermi liquid limit at low enough temperature and 
sufficiently large antinesting or low superconducting Tc. 

Quantum critical regime can be understood as an 
anomalous increase of umklapp scattering with lower- 
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ing temperature, an enhancement that turns out to be 
non uniform or anisotropic along the Fermi surface, with 
peaks at the positions where the Fermi velocity is purely 
longitudinal and nesting is relatively poor. This umk- 
lapp vertex renormalization being directly linked to an 
enhancement of SDW correlations, these can be equiv¬ 
alently considered as the source of anomalous inelastic 
scattering. The region of non Fermi liquid values of a 
are found to extend well beyond the confines of the QCP 
vicinity, a consequence of the reinforcement of umklapp 
scattering or SDW fluctuations by d-wave Cooper pair¬ 
ing. This leads to the apparent correlation of a with the 
size of Tc. 

The results have been shown to qualitatively account 
for the linear-T longitudinal resistivity observed near the 
QCP of the Bechgaard salts phase diagram and for the 
gradual growth of a toward the Fermi liquid behaviour, 
as pressure moves away from the QCP and Tc is sup¬ 
pressed. However, discrepancies between theory and ex¬ 
periment remain.The present approach to resistivity fails 
to reproduce the crossover to an apparent Fermi liquid 
in the intermediate temperature domain, namely above 
that of quantum critical effects. For systems like the 
Bechgaard salts with a simple open Fermi surface, this 
is likely to point to a mechanism of electron scattering 
that surperimposes to that of half-filling umklapp with 
increasing temperature. 
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Appendix A: The linearized Boltzmann equation 

The linearized Boltzmann equation is of the form 

E = eps ■ Vk, (Al) 

i,k' 

where the collision operator given by expression is the 
sum of four terms 




i=l 


(LNpY ^ 2 

^ k2MM 


E l\93{k±,k±2,k±3,k^4) - g3{k±,k±2,k±4,k^3)f ^ ^ ^k^k2,k3^k4^pG 


x5{el+ell-ell-eZ) 


,/0(fe2)[l-/°(fc3)][l-/°(fe4)] 


{^k,k' + ^k2,k' — ^k3,k' “ ^^^ 4 ,^')- (^ 2 ) 


In the framework of the electron gas model, we have used 

(^ 1 , ^ 21 ^ 31 ^ 3 ,^ 4 ) = 7 rhvFg 3 {k±i,k± 2 ^kps,k±Y 

for the (normalized) umklapp vertex functions evaluated 
on the Eermi surface. By separating the kinematics con¬ 
straint on the momentum conservation into its longitu¬ 
dinal and transverse componentJ^^, one has 

^k-\-k2 •,k3-\-k4-\-pG -|-/Cj_2 j^_L3+^_L4 


27r 


^_nvFS{elTellTell+el:-J:) 


(A3) 


where E = e±{k±) + ej_{k± 2 ) + ^±{k± 3 ) -h ej_{kj_ 4 ). 

The summation over the momentum vectors can be 
written as 


iNp J ‘I'kHvp Np 


LNp J 2 'Khvp Np ^ ^ 

Carrying out the integration over the energy variables 

respectively the first and the 
last two terms of ( | A2| ) , we arrive at the corresponding 
expressions entering the collision operator. 
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S lfl3(^_L, fc_L3 + fc_L4 “ ^X, ^.13, ^±4) “ 53(fcx, fcx3 + fcx4 “ ^X, fcx4, kx_^) 


P k±3,k±4 

l + e-'^^fe I]/2e^^/2 

^ 1 g/3(E/2-4) e/3E/2 

^L\' = ^J ^2 E Ifl3(fcx, fcl, fcx3, fcx4) - 53(fcx, ^^4, fcx3)r 

XiVp fc_L3,fe_L4 


X 4 j.+fe^,fcj. 3 +fe± 4 '^(^fe + ^k' ~ ^1/2) 


1 + e {el + 

1 + e^<' e^(4+<:) _ 1 ^ 




-CL\/ = -^r^ E |fl3(^x,fcx2,A:x,A:x4)-53(fc_L,fci2,fcx4,fcx)l 

Vp fcj.2,fe±4 

^ 1 + e-^-l 

^ ^k±-\-k±2,k'-\-k±4 


1 


' 1 + e-^^' + l)(e-^(^^/2-<') + 1) 


3 

E |fl3(fcx,fcx2,A:x3,fcl) -53(fcx,fcx2,fcl,fcx3)|^ 


-fc.fc' ^^2 


P /C_L2,/Cj_3 


X <5, 


1 + e-^^l 


k± -\-k±2 ,k±3-\-k'j_ 


1 + e-^<' (e/?(S3/2-4) + l)(e-^(^^/2-<') + 1) 


(A5) 


where Si, S 2 , and S 3 are obtained by substituting re- where according to (|A 1 | ) we have defined the normal- 
spectively k_\_ 2 ^ ^± 3 , and k ±4 with k'j_ in the S expression. ized quantities 7rd±0^^/vp and (p/{eEaPdj_) 

The normalized equation to solve numerically for the Carrying out the integration over e^, and rearranging the 


(p's reads 


terms we arrive at the final equation, 


i,k' 


TT^d^ T 1 


2a EfNI.> 


E 


P k'^,k^3,k^4 


2 /dS 72 

_L5 T k_\_A k_\_ , A^j_3 , k_\_A) ^35 ^P 3 T k_\_A k±_ , A^j_ 4 , ^p 3 ) I gj^q^dST^y 


2 /3S^^ /2 

T 73 (^P 5 ^P 5 ^P3 5 ^± 4 ) ^3 (^P 5 ^_L 5 ^P45 ^P3) I gj^j2h(^S^^^ 2 ^ ’^-*-^”^^-*-4 

2 jSYj'' /2 'I 

273(A:_l, A:i3,^l,^P4) - ^3(fep, fep3, fep4, 1^ gij^p(^^///2) ^fe^+fe^3,felP/c^4 (^6) 


where 


S' — e±{k±) + e_L(/c^3 + k ±4 — k_\_) + e±{k±s) + e±{k± 4 P^ 
S" = e^ik^) + e^ik'^) + e^ik^s) + e^ik^A)^ (A7) 


We shall fix the ratio d_\_/a 1.05 in ( |A 6 [ ) from the known 
values of lattice constants in the Bechgaard salts for the 
longitudinal d (7. 3A) and transverse b (7.7A) direction^. 
The solution of ( |A 6 [ ) for the 0/^^ on the Fermi surface 
leads to the evaluation of conductivity ( pT] ) as the inverse 
of resistivity. 


Appendix B: One-loop renormalization group for 
the electron gas 

1. Vertices 

The renormalization group appro ach to the quasi-ID 
electron gas is detailed in Ref^^^H^ . In the outline, the 
method starts with the segmentation of energy shells of 
thickness ^Eo{i)di^ located at ±^F^o(^) from either sides 
of the Fermi sheets, into Np patches, each centred at a 
particular value of the transverse momentum k±, where 
Eo{£) = Eoe~^ is the scaled bandwidth at step i and 
Eq = 2Ef is the initial bandwidth. The successive par¬ 
tial trace integration of electron degrees of freedom in 
the partition function leads to the renormalization or 
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flow of the scattering amplitudes gi as a function of 
At the one-loop level, the flow combines corrections from 


the electron-hole (Peierls) and electron-electron (Cooper) 
interfering channels of scattering and leads to the flow 
equations written in the compact form 


fc± 3 ,fc_L 4 )= (-251 051 +fli ° 92 + 92 ° gi)deJC.p - (51 o g 2 + g 2 0 gi)dejCc, 

deg2{k±i,k±2,k±3,k±4)= -{gi o g^ + g2 og2)deCc + {92° 92 + 93 ° 93)deCp, 

de93{k±i,k±2,k±3, k±4)= {-gi o gs - gs o gi + 92 ° 93 + 93 ° 92)deCp + 2g2 • gsde^p, (Bl) 


where = d/di. Ly^pfi are the Peierls and Cooper 
loops whose derivative at finite temperature comprises 
an integration over the patch. These take the form 




^o(^) 
4lN p 


5i 




Np dk± 
27r 


e{\Eo{£)/2 + gA,\-Eo{£)/2) 

Eo{£) + 

tanh[^£’o(^)/4] + tanh[^(£’o(^)/4 + gA^/2)] , 

(B2) 


where 


Au{k_i_,q‘]^J = - - gu£j_{guk±_ + 

+ 7?i/e_L(r?^fc_L2(4) + + e±ik±2(A)) 

(B3) 


2. Susceptibilities 


The linear response of the electron system to a source 
field of SDW and SCD order parameter leads to the cor¬ 
responding expressions of the normalized static suscepti¬ 
bility 


Xm(9m)=2 f {ffiik±)z^{kj_)fdeCf,d£ (B4) 

Je 

for fi =SDW at = (2/ci?,7r) and g =SCd at Qq = 
0.These are functions of the renormalization factor for 
the pair vertex in the corresponding channel which 
obeys the following flow equation 


dez^{k±) = {diCf,)gf,z^{k'j_), (B5) 

where 


withg^^p = k±3-k±2 = k±i-kj_4{kj_s-kj_i = k±2-k±4) 

and g_Lc = k±i -h kj _2 = ^±3 + k±4] Vp = ^ rjc = -1* 
0{x) is the Heaviside function [6>(0) = |]. 

The momentum dependence of couplings in the dis¬ 
crete convolution products 


loop variable on the right-hand side of (Bl) are 
in order g{k±,kp 4 ^k±i,k±_ - q±p)g{kp,kp 2 ^k±s,k±_ - 
q±p) for the Peierls channel, g{k±i, k_\_ 2 ^ k±, q_\_c — 
k±)g{kj_s^ k±A, k±, q±c — k±) for the Cooper channel, and 
g{kp, /c_L4, kp 2 , k± - qpp)g{kpi, k±, kps, kpi - q'^p) for 
the product of the off-diagonal Peierls channel. 


gsBW = g 2 {k'p + TT, kp, kp + TT, k'p) 

+ gsikp -7r,kp+ TT, kp, kp), 

gscd = - gi{—k'p, k'p, —kp, kp) 

- g 2 {-kp, kp, -kp, kp), (B6) 

are the combinations of momentum dependent scattering 
amplitudes that govern the flow for each susceptibility. 
The form factors assigned to the nature of ordered 
phases are /sow = 1 and fscd = cos kp . 
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